Kinetic Theory of Flocking: Derivation of Hydrodynamic Equations

It is shown how to explicitly coarse-grain the microscopic dynamics of the Vicsek model for self-propelled agents. The macroscopic transport equations are derived by means of an Enskog-type kinetic theory. Expressions for all transport coefficients at large particle speed are given. The phase transition from a disordered to a flocking state is studied numerically and analytically.Comment: 4 pages, 1 figur

Computational study of the thermal conductivity in defective carbon nanostructures

We use non-equilibrium molecular dynamics simulations to study the adverse role of defects including isotopic impurities on the thermal conductivity of carbon nanotubes, graphene and graphene nanoribbons. We find that even in structurally perfect nanotubes and graphene, isotopic impurities reduce thermal conductivity by up to one half by decreasing the phonon mean free path. An even larger thermal conductivity reduction, with the same physical origin, occurs in presence of structural defects including vacancies and edges in narrow graphene nanoribbons. Our calculations reconcile results of former studies, which differed by up to an order of magnitude, by identifying limitations of various computational approaches

Configurational temperatures and interactions in charge-stabilized colloid

We demonstrate that the configurational temperature formalism can be derived from the classical hypervirial theorem, and introduce a hierarchy of hyperconfigurational temperature definitions, which are particularly well suited for experimental studies. We then use these analytical tools to probe the electrostatic interactions in monolayers of charge-stabilized colloidal spheres confined by parallel glass surfaces. The configurational and hyperconfigurational temperatures, together with a novel thermodynamic sum rule, provide previously lacking self-consistency tests for interaction measurements based on digital video microscopy, and thereby cast new light on controversial reports of confinement-induced like-charge attractions. We further introduce a new method for measuring the pair potential directly that uses consistency of the configurational and hyperconfigurational temperatures as a set of constraints for a model-free search.Comment: 15 pages, 12 figures, submitted to J. Chem. Phy

Calculations of canonical averages from the grand canonical ensemble

Grand canonical and canonical ensembles become equivalent in the thermodynamic limit, but when the system size is finite the results obtained in the two ensembles deviate from each other. In many important cases, the canonical ensemble provides an appropriate physical description but it is often much easier to perform the calculations in the corresponding grand canonical ensemble. We present a method to compute averages in canonical ensemble based on calculations of the expectation values in grand canonical ensemble. The number of particles, which is fixed in the canonical ensemble, is not necessarily the same as the average number of particles in the grand canonical ensemble

Dynamics of a Rigid Rod in a Glassy Medium

We present simulations of the motion of a single rigid rod in a disordered static 2d-array of disk-like obstacles. The rotational, $D_{\rm R}$, and center-of-mass translational, $D_{\rm CM}$, diffusion constants are calculated for a wide range of rod length $L$ and density of obstacles $\rho$. It is found that $D_{\rm CM}$ follows the behavior predicted by kinetic theory for a hard disk with an effective radius $R(L)$. A dynamic crossover is observed in $D_{\rm R}$ for $L$ comparable to the typical distance between neighboring obstacles $d_{\rm nn}$. Using arguments from kinetic theory and reptation, we rationalize the scaling laws, dynamic exponents, and prefactors observed for $D_{\rm R}$. In analogy with the enhanced translational diffusion observed in deeply supercooled liquids, the Stokes-Einstein-Debye relation is violated for $L > 0.6d_{\rm nn}$.Comment: 8 pages, 4 figures. Major changes. To be published in Europhysics Letter

A Simple Three-Parameter Model Potential For Diatomic Systems: From Weakly and Strongly Bound Molecules to Metastable Molecular Ions

Based on a simplest molecular orbital theory of H$_{2}^{+}$, a three-parameter model potential function is proposed to describe ground-state diatomic systems with closed-shell and/or S-type valence-shell constituents over a significantly wide range of internuclear distances. More than 200 weakly and strongly bound diatomics have been studied, including neutral and singly-charged diatomics (e.g., H$_{2}$, Li$_{2}$, LiH, Cd$_{2}$, Na$_{2}^{+}$, and RbH$^{-}$), long-range bound diatomics (e.g., NaAr, CdNe, He$_{2}$, CaHe, SrHe, and BaHe), metastable molecular dications (e.g., BeH$^{++}$, AlH$^{++}$, Mg$_{2}^{++}$, and LiBa$^{++}$), and molecular trications (e.g., YHe$^{+++}$ and ScHe$^{+++}$).Comment: 5 pages, 4 figures, accepted by Physical Review Letter

Ionic Capillary Evaporation in Weakly Charged Nanopores

Using a variational field theory, we show that an electrolyte confined to a neutral cylindrical nanopore traversing a low dielectric membrane exhibits a first-order ionic liquid-vapor pseudo-phase-transition from an ionic-penetration "liquid" phase to an ionic-exclusion "vapor" phase, controlled by nanopore-modified ionic correlations and dielectric repulsion. For weakly charged nanopores, this pseudotransition survives and may shed light on the mechanism behind the rapid switching of nanopore conductivity observed in experiments.Comment: This version is accepted for publication in PR

On Calculation of Thermal Conductivity from Einstein Relation in Equilibrium MD

In equilibrium molecular dynamics, Einstein relation can be used to calculate the thermal conductivity. This method is equivalent to Green-Kubo relation and it does not require a derivation of an analytical form for the heat current. However, it is not commonly used as Green-Kubo relationship. Its wide use is hindered by the lack of a proper definition for integrated heat current (energy moment) under periodic boundary conditions. In this paper, we developed an appropriate definition for integrated heat current to calculate thermal conductivity of solids under periodic conditions. We applied this method to solid argon and silicon based systems; compared and contrasted with the Green-Kubo approach.Comment: We updated this manuscript from second version by changing the title and abstract. This paper is submitted to J. Chem. Phy

Analytical Rescaling of Polymer Dynamics from Mesoscale Simulations

We present a theoretical approach to scale the artificially fast dynamics of simulated coarse-grained polymer liquids down to its realistic value. As coarse-graining affects entropy and dissipation, two factors enter the rescaling: inclusion of intramolecular vibrational degrees of freedom, and rescaling of the friction coefficient. Because our approach is analytical, it is general and transferable. Translational and rotational diffusion of unentangled and entangled polyethylene melts, predicted from mesoscale simulations of coarse-grained polymer melts using our rescaling procedure, are in quantitative agreement with united atom simulations and with experiments.Comment: 6 pages, 2 figures, 2 table